Coexistent quantum and classical aspects of magnetization plateaux in 
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Magnetization process of ferrimagnetic Heisenberg chains of alternating spins are theoretically 
studied. The size scaling analysis with the exact diagonalization of finite systems for (S,s)=(3/2,1) 
and (2,1) indicates a multi-plateau structure in the ground-state magnetization curve for S and s 
> 1/2. The first plateau at the spontaneous magnetization can be explained by a classical origin, 
that is the Ising gap. In contrast, the second or higher one must be originated to the quantization 
of the magnetization. It is also found that all the 2s plateaux, including the classical and quantum 
ones, appear even in the isotropic case with no bond alternation. 

PACS numbers: 75.10.Jm, 75.40Mg, 75.30.Kz 



CO 



I 

C 

O 

o 



> 

o 

o 
o 



-a 
c 

o 
o 



X 



I. INTRODUCTION 



Alternating spin chains with antiferromagnetic interac- 
tions have been attracting a lot of current interest. They 
behave like a gapped antifcrromagnet based on the op- 
tical mode of the low-lying excitation, while they also 
exhibit a ferromagnetic aspect characterized by a spon- 
taneous magnetization. The coexistence of the two as- 
pects gives rise to various interesting crossover phenom- 
ena at low temperatures. (Yamamoto 2000) However, few 
quantum aspects have ever been reported for the systems. 
Then we consider an interesting phenomenon caused by 
a quantum mechanism, called quantization of magnetiza- 
tion, on the stage of the quantum ferrimagnetic chains. It 
would be observed as a plateau in the ground-state mag- 
netization curve. Recently many theoretical (Oshikawa 
et al. 1997, Totsuka 1998, Tonegawa et al. 1996, Cabra et 
al. 1997, Cabra and Grynberg 1999 and Sakai and Taka- 
hashi 1998) and experimental (Narumi et al. 1998 and 
Shiramura et al. 1998) studies suggested the realization 
of the magnetization plateaux in various systems. 

The previous works (Yamamoto and Sakai 1998 and 
Sakai and Yamamoto 1998) by the present authors indi- 
cated an important role of the quantum fluctuation to 
stabilize the plateau against the planar anisotropy in the 
ferrimagnetic chain. It results in the existence of the 
plateau even in the XY model of the mixed spins 1, and 
1/2. On the other hand, the classical mixed spin systems 
have the same plateau in the isotropic case. It implies 
that the plateau is originally based on a classical mecha- 
nism, although it is stabilized by a quantum effect. Thus 
it is difficult to say that the plateau symbolizes the quan- 
tum nature of the ferrimagnetic chains. In the present 
paper, other plateaux, essentially based on a quantum 
mechanism, are revealed to coexist with the above men- 
tioned classical plateau, when both spins are lager than 
1/2. 



II. QUANTUM AND CLASSICAL PLATEAUX 

The recent exact treatment (Oshikawa et al. 1997) for 
general quantum spin chains suggested that a magneti- 
zation plateau can appear based on the quantization of 
the magnetization under the condition, 



integer, 



(1) 



where «S U nit and m are the total spin and magnetization 
per unit cell. Note that the relation (|l|) is only a neces- 
sary condition. Thus it doesn't guarantee the existence 
of the plateau nor specify any mechanisms of its forma- 
tion. The condition (Q) is still valid for the mixed spin 
chains of S and s (S > s), described by the Hamiltonian 

N 

H = £[(1 + S)(Sj ■ 8j ) a + (1 - 6)( Sj ■ S J+1 ) a - H(S* + S |)] , 

(2) 

with (S ■ s) a = S x s x + S y sy + aS z s z . In the isotropic 
case (a — 1) with no bond alternation (6 — 0), the 
system has a spontaneous magnetization m s = S — s. 
Because of the antiferromagnetic gap the ground state 
with m s is so stable against the excitation increasing m 
that there appears a magnetization plateau at m = m s . 
(Kuramoto 1998) The previous works (Yamamoto and 
Sakai 1998 and Sakai and Yamamoto 1998) on the most 
quantized system (5, s)=(l,l/2) suggested that the quan- 
tum fluctuation stabilizes the plateau against the XY- 
like anisotropy (a < 1) and the plateau phase extends 
to the Kosterlitz-Thouless phase boundary in the ferro- 
magnetic region (a < 0). However, this plateau phase 
also includes the Ising limit (a — ► oo) without any other 
boundaries. Thus it is difficult to identify some quan- 
tum effects on the plateau formation, because it can- 
not be clearly distinguished from the Ising gap based on 
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FIG. 1. Low- lying excitation spectra calculated by the 
quantum Monte Carlo simulation (QMC), the spin wave 
theory and the perturbation from the decoupled dimer for 
(a)(3/2,l) and (b)(2,l). 

IV. VARIATIONAL APPROACH 

According to the condition of the quantization (P, the 
mixed-spin chains of (3/2,1) and (2,1) possibly have two 
plateaux at m = m s and m s + 1 . In order to characterize 
these plateaux, we introduce a variational wave function 
for the ground state of the model (||) as 

ig>= c N nf =1 (A]) 2s (6]) 2 io) 

where cn and Oy B are the mixing coefficients. Using the 
variational wave function, the ground-state phase dia- 
gram in the HS plane is obtained, as shown in Figs. 2 
(a) and (b) for (3/2,1) and (2,1), respectively, where we 
restrict us on the Heisenberg point (a = 1). On the 
first step of the magnetization process for each system, 
there exists a crossover point S c between the Neel (N) and 
double-bond dimer (DBD) states. In contrast, the second 
step toward the saturation (S) is always the single-bond 
dimer (SBD) state. These two steps before the saturation 



respectively. (Sakai and Takahashi 1998) Thus the scaled 
quantity NA N (m) is a good probe of the plateau. No 
size dependence means that the system is gapless. The 
scaled quantity of the system (3/2,1) at (a)m = 1/2 and 
(b)m = 3/2 is shown as a function of a for 5 = in Fig. 
3. Fig. 3 (a) clearly shows that opening plateau around 
a = 1 vanishes at some critical value a c for to = 1/2 and 
a gapless phase lies in the region of a < a c . The scaled 
gap in Fig. 3 (b) also indicates the existence of the sec- 
ond plateau around a = 1, although the size dependence 
is much smaller than that of the first plateau. It implies 
that the second plateau is much smaller than the first 
one. 
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FIG. 3. Scaled quantity JVZ\jv(m) versus a at m 
and m = | (b) in the case of (S, s) = (|, 1). 



(a) 



FIG. 2. The variational ground-state phase diagrams on 
the 5#-plane with a = 1 for (a)(3/2,l) and (b)(2,l). Phases 
are denoted as N: Neel, DBD: double-bond dimer, SBD: sin- 
gle-bond dimer and S: saturation, respectively. 



V. PHASE DIAGRAMS 

In order to confirm the coexistence of the two plateaux 
even for a = 1 and 5 = 0, we perform a size scaling 
analysis with the exact diagonalization of finite systems 
up to N = 12 to present the phase diagrams in the 5a 
plane. E(N, M) denotes the lowest energy in the sub- 
space with a fixed magnetization M for the Hamiltonian 
(|J) without the Zeeman term. The upper and lower 
bounds of the external field which induces the ground- 
state magnetization M are expressed as H± (TV, M) = 
±E(N,M ± 1) =F E(N,M). The length of the plateau 
with the unit-cell magnetization to = M/N is obtained 
as zljv(m) = H + (N,M) - H_{N,M). The quantity 
An(iti) also corresponds to the sum of the two excita- 
tion gaps for increasing and reducing the magnetization, 



To investigate the critical property around a c , we use 
the size scaling formula based on the conformal field the- 
ory (Cardy 1984, Blote et al. 1986 and Affleck 1986) 



1 , , s 7TW S C 

-E(N,M)~e(m)-^ 



and 



A N (m) 



N 



(5) 



(6) 



with the following notations: e(rn) is the ground state en- 
ergy of the bulk system; v s is the sound velocity derived 
from the derivative of the dispersion curve at k = 0; c 
is the central charge; r\ is the critical exponent defined 
by the spin-correlation function (s§sj:) ~ r~ n , where 
is some relevant spin operator. They are valid at gapless 
points. The calculated c and r\ of the systems (3/2,1) 
and (2,1) at m = m s and to = m s + 1 indicate the 
following properties: as a decreases from 1 with fixed 
5, the first and second plateaux vanish at the different 
Kosterlitz-Thouless critical points a c \ and a C 2, respec- 
tively, where 77 is 1/4 in common; the gapless spin fluid 
phase characterized by c = 1 lies in the region a < a c \ 
(a C 2) at to = to s (m = m s + 1). The universality of the 
phase boundary of both plateaux are the same as that of 
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the unique plateau of the system (1,1/2) (Yamamoto and 
Sakai 1999 and Sakai and Yamamoto 1999). Thus we de- 
termine the gapless-plateau phase boundary by rj = 1/4 
for both plateaux. Thus-obtained phase boundaries for 
the first and second plateaux are shown together as solid 
lines in Figs. 4 (a) and (b) for the systems (3/2,1) and 
(2,1), respectively. The phase diagrams obviously turn 
out the coexistence of the first (classical) and second 
(quantum) plateaux even at the most symmetric point 
(a = 1 and 5 = 0). They also exhibit an interesting 
feature; the quantum plateau phase is larger than the 
classical one (a c i < a c2 independently of S). It implies 
that the quantum plateau is more stable than the classi- 
cal one against the planar anisotropy. 
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1995), to convince of the existence of the second plateau. 
It is one of the most precise methods to estimate the 
Kosterlitz-Thouless phase boundary. Since the method 
detects the boundary as a level crossing point of the two 
relevant excitation gaps with the same scaling dimen- 
sion, the result does not suffer from the dominant log- 
arithmic size correction, which is quite serious for the 
Kosterlitz-Thouless transition. For the second plateau 
at m = S — s + 1 , the two relevant gaps are given by 

A =E 2 {L,M 2 )-E(L,M 2 ), (7) 

A 4 = [E(L, M 2 + 4) + E(L, M 2 - 4) - 2E(L, M 2 )]/2, 

(8) 

where E 2 (L, M) is the second eigenvalue in the same sub- 
space as E(L,M) and M 2 is defined as M 2 = (S-s + l)N. 
The two excitations have a common scaling dimension 
2. The calculated Ao and A4 of the system (3/2,1) for 
5 = are plotted versus a in Fig. 5. It suggests that the 
phase boundary is easily determined as a crossing point, 
almost independent of the system size. Thus-obtained 
boundaries for the second plateaux are shown as dashed 
curves in Figs. 4 (a) and (b). The results have a little 
deviation from the boundaries by 77 — 1/4, because the 
latter includes the logarithmic size correction. Anyway 
they lead to the same conclusion; the coexistence of the 
classical and quantum plateaux for a = 1 and 6 = 0. 



0.4 





\ two plateaux 




N, (m=l, 2) 


no plateau 


x Y J A 


(b) (2,1) 





FIG. 4. The ground-state phase diagrams on the a<5-planc 
for (a)(3/2,l) and (b)(2,l). Solid lines are the phase bound- 
aries determined by rj — 1/4. Dashed lines are the bound- 
aries of the second plateaux determined by the level spec- 
troscopy. A small difference between the two results for the 
second plateaux is due to the logarithmic size correction which 
should appear in the former method. 

Such a slight size dependence of the scaled quantity 
iVAjv(r7j) for the second plateau as shown in Fig. 3 
(b) might make us doubt its existence for a = 1 and 
S = 0. Then we perform another analysis, called level 
spectroscopy (Okamoto and Nomura 1992 and Nomura 
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FIG. 5. Ao and A 4 of the system (3/2,1) with 3 = ver- 
sus a. It indicates that the crossing point of the two gaps is 
almost independent of the system size. 

The Neel-dimcr crossover point in the first plateau in- 
dicated by the variational method in the previous section 
is not detected as any phase boundaries by these numeri- 
cal analyses. It suggests that the Neel and dimer pictures 
cannot be distinguished clearly for the first plateau. In 
fact it is trivially revealed that in the 5a phase diagram 
at m — m s the isotropic dimer point (a = 1 and 5 = 1) is 
connected to the Ising limit (a — > 00 and <5 = 0) via the 
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Ising dimer limit (a — * oo and 5 = 1) through no phase 
transition or crossover. It implies that the first plateau 
always bears an aspect of the Ising gap even for large S. 



(a) (3/2,1) 




(5,oc)=(0,0.0) 
(8,a)=(0,0.4) 
(S,<x)=(0,0.6) 
(8,a)=(0,1.0) 



H 




(S,<x)=(0,0.0) 
(S,a)=(0,0.4) 
(S,a)=(0,0.6) 
(8,a)=(0,1.0) 



FIG. 6. The ground-state magnetization curves of the 
quantum system with 8 = at various values of a for 
(a)(3/2,l) and (b)(2,l). 



only the quantum one would survive. In order to clar- 
ify the difference in the mechanism of the gap formation 
between the first and second plateaux, we also present 
the magnetization curves of the classical Heisenberg spin 
systems described by the same Hamiltonian (|J) with the 
same amplitudes (S,s) — (3/2,1) and (2,1) in Figs. 8 
(a) and (b), respectively. (The results are independent 
of 5.) The classical systems clearly have the first plateau 
in the isotropic case, while there appears no plateau cor- 
responding to the second one. It also supports the quan- 
tum nature of the second plateau. These plateaux of the 
classical systems vanishes even for a slight anisotropy; 
a c = 0.980 and 0.943 for (S, s)=(3/2,l) and (2,1), repec- 
tively. In comparison with these critical values, the phase 
boundaries in the quantum systems in Figs. 4 (a) and (b) 
imply that quantum fluctuation stabilies even the first 
plateau. Thus, in gcneal, the quantum effect is expected 
to toughen every field-induced gap against the planar 
anisotropy. Nevertheless the first and second plateaux 
should be distinguished, because the former appears even 
in the classical limit, while the latter doesn't exist until 
the spin is quantized. 



(a) (3/2,1) 




(S,oc)=(0.4,-0.2) 
(8,<x)=(0.4,0.0) 
(S,oc)=(0.4,0.6) 
(5,oc)=(0.4,1.0) 



H 



VI. MAGNETIZATION CURVE 

Finally we present the ground state magnetization 
curve in several cases of the systems (3/2,1) and (2,1). 
The curve is given by extrapolating H±(N,M) to the 
thermodynamic limit using the size scaling (Sakai and 
Takahashi 1998) based on the conformal field theory 
at gapless points and the Shanks transformation for 
plateaux. We show only the results of a suitable polyno- 
mial fitting to thus-obtained points, in Figs. 6 (6 = 0) 
and 7 (6 = 0.4), where the labels (a) and (b) indicate the 
systems (3/2,1) and (2,1), respectively. They visualize 
the coexistence of the classical and quantum plateaux 
at the Heisenberg point. These results also justify the 
above mentioned feature; the quantum plateau is smaller 
at the Heisenberg point but more stable against the XY- 
like anisotropy, rather than the cassical one. Therefore, 
if these systems lose a plateau due to the anisotropy, 




FIG. 7. The ground-state magnetization curves of the 
quantum system with 8 = 0.4 at various values of a for 
(a)(3/2,l) and (b)(2,l). 
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FIG. 8. The ground-state magnetization curves of the 
classical system at various values of a for (a) (3/2,1) and 
(b)(2,l). ' 



VII. CONCLUDING REMARKS 

The above investigations turn out the coexistence of 
the classical and quantum plateaux in the ground-state 
magnetization curve of the mixed spin chains of (3/2,1) 
and (2,1). The conclusion is easily generalized for (S,s) 
(S > s > 1/2); the magnetization curve has 2s plateaux 
and only the initial one is based on a classical mechanism, 
while the other ones originate in quantum correlations. 

In most previous works on the magnetizaion plateau 
the gap formation were based on the bond polymeriza- 
tion. (Totsuka 1998, Tonegawa et al. 1996, Cabra et al. 
1997, Cabra and Grynberg 1999) In contast, the present 
proposal of the plateau in ferrimagnetic chains is a pio- 
neering trial to explore a novel mechanism of the field- 
induced gap associated to the spin polymerization. The 
bimetallic chains such as MM'(pbaOH)(H20)3-nH20 
(Karn 1989 and Kahn et al. 1995) are good candidates 
to realize the spin plymerization. Unfortunately, most of 
them are in the case of M'=Cu, that is s = 1/2. In fact a 
few compounds with other metals were also synthesized, 
for example, MM'(EDTA)6H 2 (MM'=CoNi, MnCo 



and MnNi). However, the case of MM'=CoNi yields 
(S',s)=(l(Ni),l/2(Co)) (Drillon et al. 1985) and MnCo 
a large Ising-like anisotropy. (Drillon et al. 1986) Thus 
they are not any suitable stages to search for the quantum 
plateau. Among the series of bimetallic chains, the most 
suitable compound might be MnNi(EDTA)6H20, which 
is well described by the ID spin-alternating Heisenberg 
model of (5/2,1). (Drillon et al. 1986) The magnetization 
measurement on it would be interesting to investigate a 
possible quantum plateau at m — 5/2, as well as the 
classical one at m — 3/2. 

One of the most important remarks in the present work 
is the coexistence of the quantum and classical plateaux 
even at the most symmetric point (a = 1 and 5 = 0). 
In order to examine it, compounds consisting of metals 
and stable organic radicals (Caneschi et al. 1989 and 
Markosyan et al. 1998) might be a more ideal stage, 
because the organic radicals lead to entirely isotropic 
spin systems, rather than the bimetallic chains with in- 
evitable Ising-like anisotropy. The metal-radical complex 
also has a lot of variations. The recently synthesized one 
{Mn(hfac) 2 }3(3R) 2 (Markosyan et al. 1998) has been in- 
vestigated to realize the (5/2,3/2) spin chain. We hope 
the present calculations will stimulate not only further 
theoretical investigations, but also experimental explo- 
rations into the magnetization plateaux in ferrimagnets 
as spin-polymerized materials. 
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